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ABSTRACT 

A non-Markovian theory of population dynamics is to simulate the anti-cancer drug distribution between 
malignant and the hosting normal cell pools. The model takes into account both the cell life span and the 
proliferation rate differences. This new simulation approach looks promising for its potential to optimize a 
chemotherapeutic strategy by choosing the scheme with a higher degree of the drug-tumor selectivity. The 
pre-test designed simulation mode fits nicely the experimental data on Porphylleren-MC16 (PMC 16) phar- 
macokinetics patterns including the allometric plots revealed for this novel medicinal nanoparticle possess- 
ing some anti-cancer potential and intervening into the oxygen-independent ATP production mechanisms. 
(Int J BiomedSci 2010; 6(l):71-76) 
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INTRODUCTION 

General Background 

Drug toxicity to normal tissues and the emergence of 
drug-resistance are the major factors in limiting respon- 
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siveness to chemotherapy of cancer and many parasitic 
diseases (1-3). A non-Markovian theory of population 
dynamics in harshly varying object-surrounding environ- 
ments could be applied to chemotherapeutic systems (4, 
5). The original model presented suggests a new policy 
for improving responsiveness to the phase-specific drugs 
taking into account their non-discriminative distribution 
within a cell pool consisting of slow and fast regenerating 
(proliferating) populations. 

A Non-Markovian Population Dynamics Model to 
Propose 

The dynamics under various drug regimens of popula- 
tions that differ in life-cycle parameters is simulated using 
a computer model whose simplest form is given in eq. (a): 

x(t) = kx(t-z)[l-D(t)] (a) 
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where x(t) is the population density at time t, X is the birth 
rate, r is the generation time, and D(t) denotes the environ- 
mental process, so that D(t)>0 corresponds to the occur- 
rence of effective concentration of the drug in the system. 
Using this model the elimination time of malignant popula- 
tion (TJ and that of the limiting host population (7^) were 
estimated, and the elimination coefficient, Z, measuring the 
treatment efficacy, was calculated according to eq. (b): 
Z = 1 - Tm / Th (b) 
Simulation results (Fig. 1) show that treatment efficacy 
is a nonmonotonic function of the relation between the cell 
generation time and the period of drug administration, 
with maximal occurring when the limiting host cell cycle 
length is a multiple of the chemotherapeutic period. Ana- 
lytical results further show that in fully periodic systems 
elimination time, T, is given by eq. (c), for z>S>z/2: 




Figure 1. Simulation results of the model for two populations 
that vary in their life cycle parameters. The elimination coef- 
ficient, Z, for a malignant and a host population is calculated 
under different drug regimens represented by the effective 
period of the drug, S, and by the duration of the drug-free inter- 
val, co. Generation time of the malignant cells is 18 hr and their 
resistant life phase is normally distributed with average 6 hr 
and variance, o = 7; generation-time of the host cells is 16 hr 
and their resistant life phase is normally distributed with aver- 
age 10 hr and variance, a = 0.5; S is normally distributed with 
variance, o = 5/10, and co is constant; intrinsic birth rate is, 
X = 2; initial population size is x(0) = 5. 



T=zco/\z-(S + co)\ (c) 

In eq. (c), S is the duration of the period in which the 
drug is effective, and co is the period in which the drug 
dosage is below efficiency. The point, r = S+co, is a singu- 
lar point with T being infinite. 

To proceed the simulation data, the most common 
drug distribution and cell proliferation patterns (4, 5) were 
treated using a Sigma Plot CX600 algorithm with the La- 
bRun-06 software processed in the HP 9107-2SQ Graph 
Design Analytical Unit. 

These and previous results point out the harm in 
maintaining effective phase-specific drug dosage for 
relatively long periods: Figure 1 shows that increasing 
duration of individual drug-pulses damps down the syn- 
chronization effect, so that the ability to discriminate be- 
tween the malignant and the host populations diminishes 
greatly. Further simulations (Fig. 2) suggest that elimina- 
tion time decreases with increasing variance in cell-cycle 
length. Thus, due to larger variance in their generation- 
time (3, 6, 7), neoplastic populations are expected to be 
eliminated faster than the normal limiting populations, 
when the frequency of drug administration is a multiple 
of the average generation-time of the latter; this will be 




Figure 2. Simulation results of the model in which the average 
elimination time, T, is plotted as a function of variance in the 
generation time, r, for regimens that differ in drug's effective 
period, 5. Cell-cycle length is 12 hr; the duration of the drug- 
free interval is 7.5 hr; other parameters as in Figure 1. 
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so even if the two populations have similar average gen- 
eration time. The above results suggest that cytotoxic- 
ity to normal tissues can be reduced by a policy of short 
drug pulses whose period is a multiple of the limiting 
host population's generation time. 

That increasing variance in the temporal parameters re- 
duces the average extinction time may bear upon chemother- 
apeutic strategies of noncytotoxic drugs. Figure 2 implies 
that regimen resistance can emerge, under fully periodic 
drug regimens, by selection in the pathogen's population of 
lines whose generation-time is similar to the period of drug 
application; synchronization with the chemotherapeutic pe- 
riod offers these lines a temporal refuge from the detrimen- 
tal effect of the drug. Regimen resistance can be limited if 
artificial stochasticity is introduced in the drug regimen by 
randomizing the timing of drug administration. 

So the strategic decision making algorithm described 
is applicable to the numerous stochastic drug-cell/tissue 
distribution models (5, 6). 

Once the drug-receptor specific recognition involved 
(1, 2), this algorithm might be complicated in a way of the 
Penmann-Dalbreaux non-parametric / multi-variant dis- 
criminative approximation developments [5]. In its present 
form, nonetheless, the model proposed is surely suitable 
for all known cases where the distribution of drug recep- 
tors between malignant and the tumor- surrounding intact 
cells is about a non-discriminatory equal. 

This work is a first report to trace a direct link be- 
tween 

(a) The malignant/normal cells cytostatic drug distri- 
bution patterns, i.e. tumor-drug targeting selectivity, and 

(b) The prediction-making power of a non-Markovian 
population dynamics in its invariant version. 

This allowed to propose a simple and efficient compu- 
tational simulation method (Figures 1, 2) suitable for the 
in-advance prediction making optimization of chemother- 
apy schemes in a vast variety of cases as long as the drug 
receptor(s) distribution doesn't show a marked asymmetry 
in a mixed pool of differentiated and malignant cells. An 
obvious advantage of the model proposed is that it takes 
into account an easily available data on the cell life span 
values and proliferation rates. 

Porphylleren - MC16 or, in brief, PMC 16 (Figure 3), a 
recently developed pharmacophore based on the ferropor- 
phyrin-adducted cyclohexyl-fullerene - C 60 "ball" (8, 9), 
is now proven of being safe and efficient nanodrug capable 
to correct some oxygen independent substrate phosphor- 
ylation disorders and also to modulate the "chaotically 
fluctuating" ATP / GTP production processes in aggres- 



sively proliferating mammalian cells (10, 11). This novel 
nanodrug, therefore, possesses some anti-cancer activity 
potential (12-14). 

So PMC16 could be treated as a subject for the non- 
Markovian pharmacokinetics model (Figures 1, 2). This 
might be further developed towards an extrapolation of 
animal data to predict pharmacokinetic parameters by 
allometric scaling which is an often-used tool in drug 
development with multiple approaches available at vari- 
able success rates. In the most frequently used approach, 
pharmacokinetic parameters between different species are 
related via body weight using a power function: 

P = a • W b (d) 
where P is the pharmacokinetic parameter scaled, Wis the 
body weight in kilograms, a is the allometric coefficient, 
and b is the allometric exponent, a and b are specific con- 
stants for each parameter of a compound. General tenden- 
cies for the allometric exponent we estimated for PMC16 
are 0.75 for rate constants (i.e. clearance, elimination rate 
constant), 1 for volumes of distribution, and 0.25 for half- 
lives (see Results & Discussion, Fig. 6). 

For most traditional, small-molecule drugs, allometric 
scaling is often imprecise, especially if hepatic metabo- 
lism is a major elimination path way and/or if there are 
interspecies differences in metabolism (13). However, all 
hardly metabolized xenobiotics with molecular mass val- 



coo- -ooc 




Figure 3. Porphylleren - MC16 (PMC16) (Buckminsterfuller- 
ene (C 60 ) - 2 - (butadiene - 1 - yl) - tetra (0 - y - aminobutyryl 
- o - phthalyl) -4-cyclohexyl-2-ferroporphyrin). 
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ues higher than 2.0 kDa and, particularly, the metaboli- 
cally stable carbon nanoparticles and their water-soluble 
derivatives were found to be perfectly suitable for allome- 
tric scaling "on the edge" between pre-clinical and clinical 
trial steps (9, 11, 12, 14-16). 

METHODS 

For the in vivo administration, whatever animal (mouse, 
rat, rabbit, dog, goat) employed, a single intravenous injec- 
tion of PMC16 dissolved in sterile physiological solution 
has been made as described in (9, 11). All animals were 
kept on a standard vitamins enriched diet and were starv- 
ing for 24 hrs prior to experiment. 

To estimate all major conventional pharmacokinetics 
patterns, a routine procedure specified in (11, 12) has been 
employed. In brief, the PMC16 identification/quantification 
in any biomaterial tested (urine, blood cells and plasma, 
tissue homogenates) was carried out using the lyophilized 
CS 2 /acetone-extracts treated by an original Trident-SLC4 
/ HPLC - MALDI chromato-mass-spectrometry technique 
in the Varian SQ400 LC-MS Analyser (10, 11). 

The reagents used, all of Analytical grade, were pur- 
chased from the Bio-Rad Corp. Moscow Division, Mos- 
cow, Russia. 

A non-parametric statistical treatment (for "n" equal 
to 6 or lower) was performed using the Sigma Biostat A6 
software package to evaluate the significance of experi- 
ment/control differences. A Penmann-Dalbreaux approxi- 
mation technique has then been employed once the stan- 
dard error of the mean (SEM) values were found to be not 
higher than 6.5% of the mean (0.020-0.065 M limit range) 
whatever pair compared. To manage this, the HP LabRun- 
06 software operating graph-design algorithm has been 
processed in the HP6100-J2A Analytical Unit to reach the 
perfection for all graphics build. 

RESULTS AND DISCUSSION 

Pharmacokinetics / pharmacodynamics simultaneous 
(PK/PD) modeling does not only allow for a continuous de- 
scription of the time course of effect as a function of the dos- 
ing regimen and comprehensive summary of available data 
but also enables testing of competing hypotheses regarding 
processes altered by the drug, allows to make predictions of 
drug effects under new conditions, and facilitates the esti- 
mation of inaccessible system variables (Figure 4). 

The application of PK/PD modeling is beneficial in all 
phases of preclinical and clinical drug development, with 



a focus on dosage optimization and identification of co- 
variates that are causal for intra- and interindividual dif- 
ferences in drug response and/or toxicity. It has recently 
further been endorsed by the publication of the Exposure- 
Response Guidance document by the U.S. Food and Drug 
Administration (14). Mechanism-based PK/PD modeling 
appreciating the physiological events involved in the elab- 
oration of the observed effect has been promoted as supe- 
rior modeling approach compared to empirical modeling, 
especially because it does not only describe observations 
but also offers some insight into the underlying biological 
processes involved and thus provides flexibility in extrap- 
olating the model to other clinical situations (15). 

While drug concentrations are usually analytically 
quantified in plasma, serum, or blood, the magnitude of 
the observed response is determined by the concentration 
of the drug at its effect site, the site of action in the target 
tissue. The relationship between the drug concentration 
in plasma and at the effect site may either be constant or 
undergo time -dependent changes. If equilibrium between 
both concentrations is rapidly achieved or the site of ac- 
tion is within plasma, serum or blood, there is practically 
a constant relationship between both concentrations with 
no temporal delay between plasma and effect site. In this 
case, measured concentrations can directly serve as input 
for a pharmacodynamic model. The most frequently used 
direct link pharmacodynamic model is the sigmoid E max - 
Model: 
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Figure 4. Pharmacokinetics/Pharmacodynamics (PK/PD) 
Combined analysis general scheme. 
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With E as minimum achievable effect, C as drug 

max ° 

concentration at the effect site, and EC 50 the concentra- 
tion of the drug that produces half of the maximum effect. 
The Hill-coefficient n is a shape factor that allows for an 
improved fit of the relationship to the observed data. Thus, 
a direct link model directly connects measured concentra- 
tion to the observed effect without any temporal delay. As 
seen from Figure 5, this was observed indeed in PMC16 
pharmacokinetics experiments while the dose-interval re- 
gime parameters were in fact chosen for a drug-lympho- 
cyte discrimination task pre-solved due to a non-Markov- 
ian algorithm employed earlier (Figures 1, 2). 

All allometric plots obtained for the key PMC16 phar- 
macokinetics parameters (Figure 6) reveals the regular- 
ity surely positive for a clinical trial planning to come 
(16, 17). 

A two-compartment model proposed (PMC16) is fitted 
to the following non-Markovian compatible pharmacoki- 
netics data with both inter-individual and inter-occassional 




Figure 5. Nonlinear pharmacokinetics of pmcl6, presented 
as measured and modeled plasma concentration - time curves 
(Mean ± SE) after a single intravenous injection of 0.1 mg kg 1 
(n = 5), 1.0 mg kg 1 (n = 3), and 10 mg kg 1 (n = 8) in rats. 




Figure 6. Allometric plots of the pharmacokinetic parameters: clearance, volume of the central compartment (v), volume of distribution 
at steady state (vj, and elimination half-life of pmcl6. Each data point within the plot represents an averaged value of the pharmacoki- 
netic parameter with increasing weight from mouse, rat, rabbit (3.5 kg), dog (6.3 kg), and goat, respectively. The solid line is the best fit 
with a power function to relate pharmacokinetic parameters to body weight. 
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random effects on CL, V, Q, and V2 corrected to a error 
best described the pattern of residual error (17-19). 
Inter-individual variability with covariate model: 

A. Pharmacokinetic model 

D r(a-£) (gx0+ ()i-g ex P -(fixt)\(f) 

B. Covariate model 

CL. = [0 3 x OCC 1 + 6 4 (HT-75)] exp(r, a ) (g) 
V . = [0 X - (GFR -80) x 0 2 ] x exp(^) (h) 
k 12j = (^)xexp(^) (i) 

C. Population pharmacokinetic model (j) 

D. Population parameters 
0 = 15.5 

0= 0.229 

Thus, this peculiar segment of our present study is all 
about to get the pre-clinical trial (animal tests) derived in- 
formation on the drug blood plasma level in humans which 
would be, in turn, minimally enough to provide an optimal 
drug-cell targeting distribution we may expect (predict) 
following our non-Markovian path in a chemotherapy 
strategy developments. 

So the simulation approach we proposed would improve 
the efficacy of some new anti-cancer drugs pharmacoki- 
netics/pharmacodynamics study at the edge of pre-clinical 
and clinical trials. This approach provides an opportunity 
to save time and resources owing to a simple prediction- 
power model that reveals the drug distribution among the 
targeted compartments with different turnover rates. 
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